The singular limit of a haptotaxis model with 
bistable growth 



Elisabeth LOGAK * and Chao WANG 



Abstract 

We consider a model for haptotaxis with bistable growth and study 
its singular limit. This yiels an interface motion where the normal 
velocity of the interface depends on the mean curvature and on some 
nonlocal haptotaxis term. We prove the result for general initial data 
after establishing a result about generation of interface in a small time. 



1 Introduction 

In this article we consider a model for haptotaxis with growth. Haptotaxis 
is the directed motion of cells by migration up a gradient of cellular adhe- 
sion sites located in the extracellular matrix (ECM). This process appears 
in tumor invasion and is involved in the first stage of proliferation. It also 
plays an important role in wound healing. 

The basic mechanism involves 3 main cellular components: the tumor cells, 
the Extracellular Matrix (ECM), and some Matrix Degrading Enzymes 
(MDE). Tumor cells migrate in response to gradients of some ECM pro- 
teins. Those ECM proteins are degraded by MDE, those enzymes being 
produced by tumor cells themselves. Moreover, both tumor cells and MDE 
diffuse in the cellular medium but ECM proteins do not diffuse. 
This mechanism is reminiscent of chemotaxis, which is accounting for the 
directed migration of biological individuals (e.g. bacteria) towards higher 
gradients of some chemical substance. Chemotaxis often works 
gregating mechanism, which is reflected in the blow-up of solutions of the 
Keller-Segel model, a phenomenon that has been widely studied in the re- 
cent years. However there is a major difference between chemotaxis and 
haptotaxis: since ECM proteins do not diffuse, instead of the elliptic or 
parabolic coupling appearing in chemotaxis, the haptotaxis model involves 
an ODE coupling between the concentration of ECM proteins and the MDE 
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concentration. This is also the case in angiogenesis model (cf (8]) but mod- 
els for haptotaxis involve (at least) 3 equations, whereas angiogenesis is a 
coupled system of 2 equations. 



We now give a brief review of the mathematical literature related to hap- 
totaxis modelling. The relevant variables are the tumor cells concentration, 
the Extracellular Matrix concentration (ECM) the Matrix Degrading En- 
zymes concentration (MDE) as well as the oxygen concentration. A hybrid 
model using PDEs and cellular automata has been proposed by Anderson 
[3], involving 4 components: tumor cells, ECM, MDE + Oxygen. Global 
existence for Anderson's model has been established in |14j (Walker, Webb 
(07)). Our model is a simpler version from this model involving 3 compo- 
nents where we introduce a bistable nonlinearity to model the role of changes 
in oxygen concentration. A similar model of haptotaxis with a logistic non- 
linearity is studied in |10| (and the references therein) and global existence 
is proved. Finally Chaplain, Lolas (see [6] and the references therein) pro- 
posed a combined chemotaxis-haptotaxis model with logistic source. Recent 
results by Y. Tao, M. Wang [llj and Y. Tao [12] show global existence for 
this model in dimension N < 2. Complex patterns in haptotaxis models are 
obtained numerically in |14j . and also in [15) and [3]. 

Our starting point is the haptotaxis model proposed in [T3]. In this paper, 
the authors prove global well-posedness for a large class of initial data, a 
result which strongly emphasizes the difference with Keller-Segel chemotaxis 
model. Here we consider a different version of this model, where we do 
not explicitely consider the oxygen concentration as a variable. Instead we 
replace it by a bistable nonlinearity in the equation for the cell concentration. 
Next we show that in the limit e — > 0, the solutions converge to the solutions 
of a free boundary problem where the interface motion is driven by mean 
curvature plus an haptotaxis term. 
More precisely, we study the initial value Problem (P e ) 



where is a smooth bounded domain in W N (N > 2), @t = fix [0,T] with 
T > 0, v is the exterior normal vector on 80 and A > 0, a > are strictly 
positive constants. 



u t = Au-X7-{uV X {v)) + -f{u) 



in x (0, T] 

in x (0, T] 
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The haptotaxis sensitivity function x 1S smooth and satisfies 

Mv > 0, X (v) > 0, x» > 0. 
The growth term / is bistable and is given by 

VuGM, f(u) =u(l-u)(u- -) 

so that Jq f(u)du = 0. 

We make the following assumptions about the initial data. 

1. uq, vq and mo are nonnegative C 2 functions in and we fix a constant 
Co > 1 such that 

|K|| C 2(0) + IMIc2(0) + ll m o||,72(0) < C . (1.1) 

2. vq satisfies the homogeneous Neumann boundary condition 

^ = on 80. (1.2) 
av 

3. The open set 0o defined by 

00 := {x 6 0,n o (x) > 1/2} 
is connected and 0o CC 0. 

4. To := 90o is a smooth hypersurface without boundary. 

With these assumptions 0o is a domain enclosed by the initial interface To 
and 

uo > 1/2 in O , < uo < 1/2 in \ 0„. 

The existence of a unique nonnegative solution (u £ , v £ ,m £ ) to Problem (P £ ) 
is established in Section 2. Note that it follows from (11.21) that 



T?- = Oond0x[O,r]. (1.3) 

We are interested in the asymptotic behavior of {u £ ,v £ ,m £ ) as e — >• 0. The 
asymptotic limit of Problem (P £ ) as e — > is given by the following free 
boundary Problem (P°) 



(P° 



1 in t ,tj [0,T] 
in 0\0 t ,t € [0,T] 
w ° = -AmV in x (0, T] 



m°(M) = X0t(») 



- aAm° + u° -m° in x (0, T] 



t ! 

V n = -(N-1)K + ^^ onT t = d0 t ,te(O,T] 

on 

r*]t=o = To 
u°(x,0) = u (ic) 

m°(x,0) = mo(x) x 6 $7 

8tH. 

— — = on 80 x (0,21 
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where 0t CC is a moving domain, Tt = d&t is the limit interface, n is the 
exterior normal vector on Tt V n is the normal velocity of Tt in the exterior 
direction and k is the mean curvature at each point of T t - We first establish 
the well-posedness of Problem (P°) locally in time in Section 3. Our main 
result is to prove rigorously the convergence of (u £ ,v £ ,m £ ) to (u ,v°,m ) 
for initial data satisfying the above assumptions. In a first step, we establish 
the following generation of interface property. 

Theorem 1.1 Assume that (uo,Vo,mo) satisfy the hypotheses 1-2-3-4. Let 

< rj < 1/4 be an arbitrary constant and define fj, = f (1/2) = 1/4. Then 
there exist Eq > and Mq > such that, for all e G (0, £q\ and all t G [t*,T] 
where t* = /i _1 e 2 | lne|, 

(a) for all x £ 0, we have 

< u 6 (x,t) < 1 + r?; 

(b) for all x G such that \uq{x) — \ \> Mqe, we have 

if uo(x) > — + Mqe, then u £ (x,t) > 1 — rj, 

ifuo(x) < i — Mqe, then < u £ {x,t) < rj. 
The main result reads as follows. 

Theorem 1.2 Assume that (uo,vo,mo) satisfy the hypotheses 1-2-3-4. Let 

(u £ ,v £ ,m £ ) be the solution of Problem (P £ ) and let (v°,m , F) with T = 
(T t x {t})te[0,T] be the smooth solution of the free boundary Problem (P°) 
on [0,T]. Then, as e — > 0, the solution (u £ ,v £ ,m £ ) converges to (u°,v°,m ) 
almost everywhere in Uo<<<t((® \ x More precisely, 

limu £ (x,t) =u°(x,t) a.e. in II ((0\r t )xt), 

0<t<T 

and for all a G (0, 1), 

lun\\v £ - v°\\ cl+aA1+a)/ 2(0 T ' ) =0, 

lim \\m £ - m°|| c i +a , ( i +a)/2( T) = 0. 

We actually prove a stronger convergence result concerning u £ . 

Corollary 1.3 Assume that (uq,vo, mo) satisfy the hypotheses 1-2-3-4. Then 
for anyte(0,T], 

i- ei ja / \ f 1 / or a; G fit /-, .s 

lim u (x,t) = Y0. (x) = < _ , „ . j=r- (1.4) 

e^o v 7 aw 4 v / | o /or x G \ fit v ; 
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Moreover like in [I], we also obtain the following estimate of the distance 
between the interface Tt solution of Problem (P°) and the set 

rf :={x G 0,u £ (x,t) = 1/2}. 



Theorem 1.4 There exists C > such that 



rf C N Ce (T t ) /or < t < T, 



where H r (r t ) := {x G 0, dist(x, T t ) < r} is the tubular neighborhood ofY t 
of radius r > 0. 

Corollary 1.5 r| — )■ I\ as e — > 0, uniformly in t £ [0, T] in f/ie sense of 
the Hausdorff distance. 

The organization of the paper is as follows. In section 2 we prove some 
a priori estimates, establish a comparison principle for Problem (P £ ) and 
prove the existence of a unique global solution. In section 3 we prove the 
well-posedness of the free boundary problem (P°) and obtain the existence 
of a smooth unique solution up to some time T > 0. In section 4 we establish 
the property of generation of interface. Finally in section 5 we prove the 
convergence of the solution of Problem (P £ ) to the solution of Problem (P°). 

2 A priori estimates and comparison principle 
2.1 A priori estimates 

For a given T > and a given nonnegative function u$ G C 2 (0), we define 

I T = {«6 C°(0 T ), < it < C in T and u{x,0) = u (x)}, 

where Co > 1 is the constant defined in (jl.ip . It is convenient to rewrite 
Problem (P £ ) as an evolution equation for u with a nonlocal coefficient 
H(u) = v, namely 



where for a given function u = u(x,t) G Xt, we define H (u) = v as the first 
component of the unique solution (v,m) of the auxiliary problem 



u t = Au - V • (uVx(ff(«))) + 



in x (0,T] 
on 80 x (0,T] 



< 



u(x, 0) = uo(x), 
7T = 



(2.1) 



< 



vt = —Xmv 
rrit = aAm + u — m 
v(x,0) = v (x), 
m(x, 0) = mo(x), 



infix (0, T] 
in n x (0, T] 
x £Q 
x G fi 



(2.2) 




on 90 x (0,T]. 
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The functions vq and tuq are given and satisfy 1-2. We give below some 
a priori estimates on the solution to Problem (P £ ) and state the related 
properties of H. 

Lemma 2.1 For u £ Xt, let (v,m) be the solution of Problem (|2,2|) and let 

H : Xt — > C 2 (0t) be the operator defined by H(u) = v. Then there exists 
C > only depending on T and such that 

(a) for all (ui,U2) £ X T with < u\ < U2 in 0t, the solution (vi,mi) 
of Problem (|2.2p for i = 1, 2 satisfies 

< mi < m 2 ar ^ < V2 < Wl to 0t 
so £/«z£ i/ie operator H is nonincreasing on Xt- 

(b) for all u £ Xy, 



l m llc , l+a,(l+a)/2(0 T ) < 



CCq and sup | / Am(x, s)ds\ < CCq. 
(x,t)e0 T ^° 

(c) for all u £ £/ie function v = H{u) satisfies 

\\v\\ c o(& T ) < c o and ||Vu|| c o ( t) + ||Aw|| cO( t) < CCg. 

Proof. To prove property (a), let (-ui,^) G X T with < ui < U2 in 0^. 
Since for i = 1, 2 

(mj) t — aAm, + mj = itj > in 0^, 

with 

mi\t=o = mo > and = on 80 x (0,T], 

we deduce from the standard maximum principle that < mi < m.2 in 0y. 
Next solving the equation = —Xmv we get that 

«i(x,t) = v {x)e- x ti m ^ s)ds (2.3) 

for all (x, t) £ 0^ and i = 1,2, so that i>i > > in 0j>, which proves that 
H is nonincreasing on 

In order to prove (6) and (c), note that m satisfies the linear parabolic 
equation 

m t = aAm + u — m in x (0, T] 
m(x, 0) = mo(x) x £ (2 41 

— = on 90 x (0, 71 

8v 
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with < u < Co in &t and mo > G O,. Thus it follows from the maximum 
principle and from standard parabolic estimates that there exists a constant 
C > only depending on T and such that 

< m < Co in T , | |m| | c1+q , (1+q)/2( t) < CC . (2.5) 

In view of (|2.3p . v > and t> £ < in 0t so that 

< v(x, t) < v (x) < C for all (a?, t) G T - (2.6) 

Since for all (x, t) G 0t 

Vu(x,t) = Vwo(ic)e _A - / 'o m(a: ' s)<is - Xv(x,t)( / Vm(s,s)ds), (2.7) 

JO 

it follows that there exists C > such that 



o 



\Vv(x,t)\ <\Vv (x)\+ Xv(x,t)\ I Vm(x,s)ds\ 
Since for all (x, t) £ @t 



(2i 



<CC 2 



Av(x,t) =Au (x)e^ A io tm (^) ds -2\X7v (x).( f Vm{x,s)ds)e- X ^ m ^ d 

Jo 



r-t r-t 
+X 2 v(x,t)\ / Vm(x, s)ds\ — Xv(x,t)[ / Am(x,s)ds) 
Jo Jo 

(2.9) 

it follows that 

V(x,i)e0 T , |Au(s,t)| <CC 3 + AC | / Am(x,s)ds| (2.10) 

J o 

with C > a suitable constant. 

For any fixed x G 0, we integrate the equation mj — aAra + m = u on [0, t] 
and obtain that 

t l rt 

Am(x, s)ds = — [m(x,t) — mo(x) + / (m(x,s) — u(x,s))ds] 
o a Jo 

so that in view of ()2.5p there exists a constant C > such that 

V(x,t)G0 T , \[ Am(x,s)ds\ < CC (2.11) 
•/ o 

which completes the proof of (b). Moreover in view of (|2,10p and (|2,lip . we 
conclude that there exists C > such that 

V(x,t)G0 T , \Av{x,t)\ < CCl 

and obtain the property (c) , which completes the proof of Lemma 12.11 



2.2 Existence of a global solution to Problem (P £ ) 

We prove the existence of a unique solution {u £ ,v £ ,m £ ) to Problem (P £ ) on 
0T for e > small enough. 

Lemma 2.2 Assume that (uo,vo,mo) satisfy the hypotheses 1-2-3-4. Then 
there exists £q > such that for all < e < Eq, Problem (P £ ) has a unique 
solution (u £ ,v £ ,m £ ) on x [0, T] for any T > 0. This solution satisfies 
< u £ < C in T - 

The above lemma is similar to Lemma 4.2 in [3] and we just sketch the proof. 
It relies on Schauder's fixed point theorem and on the a priori estimates on 
Problem (P £ ) obtained in Lemma 12. 11 

First let T > be arbitrarily fixed and for all u 6 Xt, let v = H(u) be 
defined as above. By the estimates of v in Lemma 12.11 there exists C > 
such that 

< v < C , \Vv\ + |Au| < CCl in T . (2.12) 
Next let u be the unique solution of 

( 1 

ut = Au - V • Wx » + -o /(«) in x (0, T] 

< u(x,0) = u (ac) xGfi (2.13) 

du „ _ 

— = on 50. 

The key point of the proof is to show that for < e < £o small enough, we 
have 

< u < C in T . 

This follows from the fact that Co is a supersolution for equation (|2.13p for 
e > small enough. Precisely, using that /(Co) < since Co > 1 and (|2.12|) . 
we have that 

C„A(x(u)) - 1/(C7 ) 
= CoCV^ACt;) + x"(«)|V(t;)| 2 ) - ^/(C ) 
> -2C 4 - 1/(C ) > 

for e > small enough. Moreover u £ C a,a//2 (0r) for some a G (0,1). 
Hence it — )• {t maps into itself and defines a compact operator. A fixed 
point of this operator obtained by Schauder's theorem is then a solution to 
Problem (P £ ). The uniqueness of solution follows from the a priori estimates 
on Problem (P £ ). For the details of the proof, we refer to [1] and [9]. 
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2.3 A comparison principle for Problem (P £ ) 

We first recall the definition of a pair of sub- and super-solutions similar to 
the one proposed in [3J. 

Definition 2.3 Let (u~ ,uf) be two smooth functions with < u~ < uf in 

Qt and q 6 < q 6 on 80 x (0,T). By definition, (u~,uf) is a pair of 
sub- and super-solutions in £1? if for any v = H(u), with u~ < u < uf in 
0T, we have 

L v [u~\ < < L v [u+\ in0 T , 
where the operator L v is defined by 

L v [4>] = <t>t ~ + V • (<PV X (v)) - \f (</>). 

Note that in Lemma 2.2, (0,Co) is a pair of sub- and super-solutions of 
Problem {P e ). It is then proved in [4] that the following comparison principle 
holds. 

Proposition 2.4 Let a pair of sub- and super- solutions (u~,uf) in 0t be 
given. Assume that 

Vx G 0, u~(x,0) < u (x) < uf(x,0), 

with (uo,vo,mo) satisfying the hypotheses 1-2. Then there exists a unique 
solution (u £ , v £ , m £ ) of Problem (P £ ) with 

V(x,t) G 0t, u7(x,t) < u £ (x,t) < uf(x,t). 



3 Well-posedness of Problem (P°) 

We establish here the existence and uniqueness of a smooth solution to the 
free boundary Problem (P°) locally in time. 

Theorem 3.1 Let Tq = 80q, where 0q CC is a C 2+a domain with 
a £ (0,1). Then there exists a time T > such that Problem (P°) has a 
unique solution (y° ,m° ,T) on [0, T] with 

T=(T t x {i}) te[0 ,T] G C 2+a ^' 2 and v\ G C 1+a ^' 2 

This theorem is similar to Theorem 2.1 in [3] and is using a contraction 
fixed-point argument in suitable Holder spaces (see Section 2 in [3]). We 
show here how it can actually be obtained using the result established in 
Theorem 2.1 in [3] and some additional properties that we state and prove 
below. 

First we introduce some notations as in [3] . We assume that Tq is parametrized 
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by some smooth (N — l)-dimensional compact manifold M without bound- 
aries which divides M. N into two pieces. We denote by N(s) the outward 
normal vector to M at s G M and define 

X:Mx(-L,+L) -> R N 

(s,s N ) h-> X(s,sat) 

where 

X(s, Sat) = s + snN(s). 

If L > is chosen small enough, X is a C^-diffeomorphism from x 
(— L, +L) onto a tubular neighborhood of .M that we denote by M. L . We 

L_ 

assume that r\) C M 2 and is given by 

T = {X(s, s N ), s N = A (s), seM} 
and that fio is the connected component of Q \ To which contains 

{x = X(s,sn), sn < A (s), s G M}. 
According to the regularity hypothesis on Tq in Theorem 13.11 Ao is a 

C 2+a 

function with 

MA II L 

\\ A o\\c°(M) < 2' 

Let T > be a fixed constant that will be chosen later. We parametrize the 
interface T = (Tt)t^[o,T] as follows 

T t = {X(s, s N ), s N = A(s, t), s G M}, (3.1) 

where A : Ai x [0,T] — > (—L,+L) is a function. By definition, we will say 
that T is C mJrCX,w ^ £L if the function A satisfies 

A G C m+a ^{M x [0,T]) 

For any function v(x,t) defined in 0t, we consider the restriction of v and 
of Vv on the interface T and we associate to v the functions w(s,t) and 
fr(s, £) defined onXx [0, T] by 

u/(s, t) = v(X(s, A(s, t)), t), (3.2) 
h(s, t) = Vv(X(s, A(s, t)),t). (3.3) 

Next we split Problem (P°) into two subproblems (p a ) and (pb), where 
Problem (p a ) is given by 

(p a ){ Vn = ~( N '-l)K + x'Wh-ftcm T t = 80 t , t€ (0,T] ^ 
1 rd t= o = To 
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and Problem is given by 



(Pb) < 



)® = —Xm°v° 

aAm° + m° 

) =0 

ou 

u°(x,t) = X0 t (x) 



dmP 



in x (0, T] 
in x (0, T] 

on 80 x (0,T] 



1 in 0t,t_G [0,7] 
in 0\0 t ,t€ [0,T] 



(3.5) 

Note that the difference between the free boundary problem in [lj and here 
concerns Problem (fib)- Let us consider 



\/(x,t)eQ T , M(x,t)= / m°(x,s)ds 



(3.6) 



The restrictions of M and VM on T are denoted a(s,t) and b(s,t) and 
defined on .M x [0, T] by 

a(a, t) = M{X(s, A(s, t)),t), (3.7) 
b(s, t) = X7M(X(s, A(s, t)), t). (3.8) 

Note that using (|2.3p and (|2.7|) we have that 

w( S ,t)=v (X(s,A(s,t)))e- Xa ^ 

and 

h(s,t) = Vv (X(s,A(s,t)))e- Xa ^ - Xw(s,t)b{s,t), 

so that w has the same regularity as a and h has the same regularity as b. 
We deduce from Problem (pt>) that M satisfies 



—aAM + M = 5(ac,t) 



in x (0, T] 
on d& x (0,T], 



(3.9) 



where 



g(x,t) = / tt°(x, s)ds + mo(x) — m°(x, t). 



The same problem 13.91 has been considered in [5j but with a right-hand-side 
g = u°. Here the function g(x,t) is continuous in time, its regularity being 
the one of a time-integral of u°. Thus we can use Theorem 2.2 in [3] and 
obtain (at least) the same regularity for (a, b) in the case considered here. 

Lemma 3.2 Let T = (T t x {t}) t e[o,T] be given by (|3.ip with 

A e C m+a '^{M x [0,T]) 
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for some m G N, m > 2 and a G (0, 1). Let M satisfy f|3.9|> and Zei a and 6 
6e associated to M by flfl.Tp and ( E£ ffj) respectively. Then 

a e c-m+a,,^^ x [0)T]) 

and 

6 G [r^'^tM x [0,T])] n /or a// < a' < a. 

By the argument in [3] we know then that Problem (p a ) defines a mapping 
(w, h) — > A and Problem (pj) defines a mapping A — )• (k;, /i) with the proper 
regularity in Holder spaces. Therefore the composition of these two map- 
pings defines a contraction in some closed ball for T > small enough. The 
unique fixed point of this contraction is the solution to Problem (P°) on 
[0,T]. This completes the proof of Theorem 13.11 

4 Generation of interface 

In this section we establish the rapid formation of transition layers in a 
neighborhood of To within a very short time interval of order e 2 | lne|. The 
width of the transition layer around To is of order e. After a short time the 
solution u £ becomes close to 1 or except in a small neighborhood of To- It 
reads precisely as follows. 

Theorem 4.1 Let uq satisfy the assumptions 1-2-3-4. Let < rj < 1/4 and 
define fi = f (1/2) = 1/4. Then there exist Eq > and Mq > such that, 
for all e G (0, eq] and t* = [i~ 1 £ 2 \ hie|, 

(a) for all x G 0, we have 

-i] < u £ (x,t*) < 1 + r); 

(b) for all x G such that \uq{x) — \\> Mqe, we have 

if uq(x) > - + Mqe, then u £ (x, £*) > 1 — rj, 

if uq(x) < — — Mqe, then u 6 (x,t*) < n. 

The above theorem relies on the construction of a suitable pair of sub- and 
super-solutions involving the solution of the bistable ODE. We refer to the 
proof of Theorem 3.1 in [1] in the simple case 5 = 0. 

5 Convergence 

We split the present section into 2 parts. In a first step we establish the 
convergence of the u £ to u° and prove Corollary ll.3p . In a second step we 
prove Theorem 11.21 as well as Theorem ll.il Theorem 11.41 and Corollary 11.51 
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In what follows, we construct a pair of sub- and super-solution uf for Prob- 
lem (P £ ) in order to control the function u £ on [t*,T]. By the compari- 
son principle it then follows that, if u~(x,0) < u £ (x,t*) < uf(x,0), then 
u7{x,t) < u £ (x,t + t*) < uf(x,t) for all (x,t) G 0t- As a result, if 
both uf and u~ converge to u°, the solution u £ also converge to u° for 
all (x,t) G T \T. 

5.1 Construction of sub- and super-solutions 

Before the construction, we present the definition of the modified signed 
distance function which is essential for our construction of sub- and super- 
solutions. Let us first define the signed distance function. 

Definition 5.1 Let T = Uo<t<T(^* x *) ^ e ^ e solution of the limit geometric 
motion Problem (P )- The signed distance function d(x,t) is defined by 

d(x,t) = { *? {x ' f T i , f° rx ^\^ (5.1) 
y J { -dist(x,T t ) forxEQt, 

where dist(x,T t ) is the distance from x to the hyperface T t in 0. 

Note that d(x,t) = on T and that \Vd(x,t)\ = 1 in a neighborhood ofT. 

In fact, rather than working with the above signed distance function d(x, t), 
we need a modified signed distance function d defined as follows. 

Definition 5.2 Let do > small enough such that d(x,t) is smooth in 

{(x,t) G 0x [0,T},\d(x,t)\ <3d } 
and such that for all t G [0, T], 

dist{T t ,d&) > 4d . 
We define the modified signed distance function d(x, t) by 

d(x,t) =((d(x,t)), 
where ((s) is a smooth increasing function on M. N defined by 

{s if \s\ < 2do 

-3d if s < -M (5.2) 
3do if s > 3do- 

Note that |Vd| = 1 in the region {\d(x,t)\ < 2d ,(x,t) G x [0,T]}. It 
follows that at x G T t , the exterior normal vector is n(x,t) = Vd(x,t), 
the normal velocity is V n (x,t) = —df(x,t) and the mean curvature is K = 
j^Ad(x,t). Therefore the motion law on T* given by Problem (P°) reads 

d t - Ad + Vd.V X (v°) = on r t = {x G fi d(x,t)=0}. (5.3) 
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By Theorem 13.14 the interface I\ is of class C 2+a ' a and v° is of class 
(ji+a , + £ Qr an y qj^q, g (g, 1), all the functions dt, Ad, Vd are Lipschitz 
continuous near I\ and Vx(v ) is continuous near IV Therefore from the 
mean value theorem applied separately on both sides of I\, it follows that 
there exists Nq > such that 

V(x,i)efi T , \dt ~ Ad + X7d.X7 x {v°)\ < N \d(x,t)\. (5.4) 

Note also that by construction, Vd(x,t) = in a neighborhood of 80. 
As in PP, the sub- and super-solutions uf are defined by 

ut = Uo( dix ' t) * £pit) )±l(t), (5-5) 

where J7o(^) is the unique solution of the stationary problem 

K + f(U ) = , , 

Uo(-oo) = 1, Uq(0) = \, *7 (+oc) = 1 ' ; 

and 

p(t) = - e -W £ ' 2 +e Lt + K 
q(t) = a(f3e-W £2 +e 2 Le Lt ) 
with L > and -fT > 1 to be chosen later. 

First note that q = e 2 o~pt, then remark that for Problem (|5.6p the unique 
solution Uq has the following properties. 

Lemma 5.3 There exist the positive constants C and A such that the fol- 
lowing estimates hold: 

< U (z) < Ce~ A|z| for z>0, 

< 1 - U (z) < Ce" A|z| forz<0. 

In addition, Uq is strictly decreasing and \Uq{z)\ + \Uq (z)\ < Ce~ A ' z /or all 
zGi 

The proof of Lemma 5.4 is given in [3]. We also note that 
u-(x,t) <U (^^-) <ut(x,t) 

and that p(t) is bounded for all < e < Eq and t E [0,T], lim £ _^o l{t) = 
for all t > 0. Therefore it follows from the definition of uf(x,t) that for all 
t€(0,T], 

limu±(M) = X0t (x) = {l j£ ^ (^j g \ Q t ( 5 - 7 ) 
The key result of this section is the following lemma. 

Lemma 5.4 There exist f3 > 0, a > smc/j £/iai /or a// > 1, we can find 
eo > and L > suc/i that for any e £ (0,£o); ( u e > uf ) is a pair of sub- 
and super-solutions for Problem (P £ ) in x [0, T]. 
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5.2 Proof of Lemma 15.41 

First note that for all (x, t) G 0t, 

u-(x,t) < U (^-) - q(t) < Uo (^-)+q(t) < u+(x,t). 

8u^~ 

Next since Vd = in a neighborhood of 80, we have that £ = on 

av 

80 x [0, T]. Let d be such that v = H(u) with u~ < u < in 0t, we show 
below that 

L v [u-} < < L„[u+], 
where the operator L v is defined by 

LM = <p t - A4> + V(0V X W) - 

Here we just consider the inequality L„[ii+] > 0, because the proof of the 
other inequality L„[u~] < is obtained by similar arguments. A direct 
computation gives us the following terms 

(uf) t = U' {^ - p t ) + q t , 

v< = u —, 

A + rr"\^d\ 2 ,Ad 

where the value of the function Uq and its derivatives are taken at the point 
d(x, t) — spit) 

- . Moreover the bistable function has the expansions 

£ 

f{ut) = f{U )+qf\UQ)+ l - q 2 f"{e), 

where 6(x,t) is a function satisfying Uq < 9 < uf. Hence, combining all the 
above, we obtain that 

L v [u+] = {ut)t-Au++Vu+Vx(v)+u+Ax(v)--2f(u+) = E 1 +E 2 +E 3 +E 4 
where 

Ei = ~q[f(U ) + \qf"{0)] - U' oPt + q t , 



U, 



E 2 = -§-(1 - \Vd\ 



2-< 



E 3 = ^(d t -Ad + Vd-V X (vo)), 
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£ 4 = _o V d • V( X (v) - x(v°)) + u+Ax(v). 

We first need to present some useful inequalities before estimating the four 
terms above, this step is exactly the same as in pQ. 

Since / (0) = /'(l) = — - , we can find < b < 1/2 and m > such that 

if U (z) G [0, b] U [1 - 6, 1] then f'{U (z)) < -m. 

Furthermore, since the region {z G W,Uq(z) G [b, 1 — b]} is compact and 
[/q < on R, there exists a constant ai > such that 

if £7o0z) G [6, 1 - 6] then U' (z) < -a l . 

Now we define 

F= sup (|/(z)| + |/(*)| + |/" (ar)!), 

-1<2<2 



/5 = j, (5-8) 

and choose a which satisfies 

< a < min(cxo, <J\, u-z)-, (5.9) 

a\ 1 4/3 , . , 

where a"o = > fi = 7^ j o"2 = — 77^ r- Hence we obtain that 

m + F p + 1 F(p + 1) 

Vz G M, -U' (z) - af'(U (z)) > Aa/3. 

Now we have already chosen the appropriate (3 and a. Let K > 1 be 
arbitrary, next we prove that L v e [uf] > provided that the constants £o > 
and L > are appropriately chosen. From now on, we suppose that the 
following inequality is satisfied 

e 2 Le LT < 1. (5.10) 

Then given any e G (0, £o)> since < a < cri, we have < q(t) < 1 for all 
t > 0. Since < U < 1, it follows that for all (x,t) G 0r 

-l<«f(x,t)<2. (5.11) 

We begin to estimate the four terms Ei, E2, E3 and E4. The estimates of 
the terms E\, E2 and E3 are similar to the estimates in [1] and we obtain 
that 

El > ^ e ~^ 2 + 2<T/3Le Lt = + d x Le L \ 
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where d = aft 2 , C x = 2a are positive constants. 

1 f\C 

\E 2 \<——^(l + \\Vd\\l) = C 2 , 
(eAdo) 

where C and A are the constants that we choose in Lemma 5.4, so that C 2 
is also a positive constant. 

We remark that in the estimate for E 2 in [I], the following assumption holds: 

e ^ + K<A. (5.12) 

For E%, we use (|5.4[) and obtain that 

\E 3 \ < C 3 (e Lt + K) + C' 3 , 

where C3 = NqC and C 3 = — - — with C and A the constants given by 

A 

Lemma 5.4 . 

Then we consider the term £4. We should know the estimates of V(x(u) — 
x{ v °)) an d Ax(v). In fact, for this term, we have the following lemma. 

Lemma 5.5 Let u be any function satisfying 

u~ < u < uf in @t 

and let (v,m) be the corresponding solution of Problem (|2.2|) with v = H{u). 
Then there exists C > depending on T and such that for all (2, t) S 0t, 

\v(x, t)\ + \Vv(x, t)\ + \Av(x, t)\<C (5.13) 

I [ (m-m°)(x,s)ds\ + \Vd(x,t) • f V(m - m°)(x, s)ds\ < Cdp%fyi) 
Jo Jo 

\(v -v°)(x,t)\ + \Vd(x,t) -V(v -v°)(x,t)\ < Cep(t) (5.15) 
where (v , m ) are given by the solution of Problem (P )- 



We prove this lemma below. Let us carry on with the proof of Lemma [57 
We write 

Vd-V(x(w)-x(« )) = x'(v)X7d-V{v-v ) + (x(v)-x{v°))Vd-X7v . (5.16) 

Since i>° is bounded in C 1+a >~2~~ for any a G (0, 1), there exists C > 0, 
such that 

ll t, °llL o °(0 T ) + II V ' i;O |Il oo (0t) - ^' 

which combined with (|5.16p . yields that 

|Vd • V( X (v) - x(v°))\ < Hx'lloolVd • V(v - v°)\ + C||Vd|| 00 ||x"||oo|u " A 

(5.17) 
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where the L°°-norms of x an d x" are considered on the interval (— C, C). 
Therefore, since x is smooth and ((Vc^oo is bounded, it follows from (|5,17p 
that for all (x, t) G 0t, there exists C > such that 

|Vd • V(x(«) - X(w°))l < Cep(t). (5.18) 

Moreover, using the smoothness of x an d the first inequality of Lemma 5.6, 
we obtain that there exists C > such that 

\Ax(v)\ < C. (5.19) 

Hence, by the above inequalities (|5.18j) . ()5.19p and the fact that \uf(x,t)\ < 
2, we obtain that for all (x, t) £ 0t, 

\E 4 \ < jCep(t) + 2C f . 

Finally substituting the expression for p and q, we obtain that there exist 
the positive constants C4, C 4 and C 4 such that 



\EA <C*+C'e-M e2 +C: 



We collect the above four estimates of E\, E2, -B3 and E4, which yield 

- C 3 (e Lt + K)-C' 3 -C±- C' 4 e~W £2 - C'ie Lt (5.20) 
Ci — e 2 C A 



£ 2 



i e -Ms' +{LC ' i _ C3 _ C >> )e Lt_ C6 



where Cq = C2 + C3K + C 3 + C4 is a positive constant. Now we set 

T 4e ' 

where Sq is small enough and satisfies the assumptions (|5.10p and (|5.12p . so 

C — e 2 C' 

that L is large enough. It also follows that 5 > and 

e z 

LC 1 — C3 — C 4 > 2^^i' 

therefore 

L v [u+] > \hC' x - C 6 > 0. 
The proof of Lemma 5.5 is now completed, with the constants /3, a given in 
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5.3 Proof of Lemma 15.51 



Lemma 15.51 gives the key estimate and is the analogue of Lemma 4.9 in [3] 
and of Lemma 2.1 in [lj. However the proof is markedly different since the 
coupling between u and v is given by a system with an ODE and a parabolic 
equation versus an elliptic equation in the two above references. 
First note that (|5.13|) is established exactly as in Lemma 12.11 (c) . 
Concerning the second inequality (|5.14p . let us recall the following properties 
of Uo given in [1]. 

Lemma 5.6 For all given a E R and z E R, we have the inequality: 

\U (z + a) - X]-oo,o](*)l < Ce~ x \ z+a \ + X ]-a,a]( Z ) 
Define w(x,t) = m(x,t) — m°(x,t), then w satisfies 



wt — aAw + w = h in @t 
^ = on80x(O,r) 



dv 

„ w(x, 0) = 0, 



(5.21) 



x 



with h = u — u° satisfying 



u u < h < ul 



e 



u° in £It- 



d(x, t) 



and 



From the definition of in ()5.5j) and from Lemma 15.61 for z 
a = ±p(t), we deduce that for all (x,t) E 0t, 

\h(x,t)\ < C(e- A l^- ! *)/-+^(*)l_ he -AMC-, t )/ s - ? ,a)| )+X{| ^ (:K t)| ^ (£)>+ ^ (t) (5-22 ) 
Let us define for all (x,t) E 0t, 

/ii(ac,t) = q(t), 



h 2 (x,t) = C (e-^y^ + e^'M^XM*,^} 



and 



M*, t) = de-W^+v (t)\ + e -xW x ,t)/e- P (t)\ )x{]d{XjmdQ} + X{|d(a . t)| < ep(t)} 

and denote by (itfi)i=i,2 3 the solutions of the three following auxiliary prob- 
lems 

(wi)t - aAwi +Wi = hi in T 

fill) ' 

(A)< ^ = ona0x(O,T) 
ov 

, Wi(x,0) =0, x E 
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Note that in view of the definition of p(t) and the inequality (j5.12|) . we have 
that for all t G [0, T] 

< K - 1 < p{t) < — (5.23) 
2e 

so that the function p is bounded away from for all t £ [0,T]. It follows 
in particular that choosing e > small enough, 

ep(t) < d /2 for all t G [0,T] 

so that \h\ < hi + h 2 + /13. Thus we deduce from the maximum principle 
that for all x £ and t G [0, T], 

\w(x, t)\ < wi(x, t) + w 2 {x, t) + w 3 (x, t). 

We now establish estimates for Wi, with i = 1, 2, 3. 
Problem (Ai) 

Set Wi(x,*) = J wx(x, s)ds, then Wi satisfies 

f (Wi)t - aAVFi + Wi = H-l in T 

^ = on90x(O,r) (5.24) 

L Wi(s,0) = 0, x G 



with, since g(t) = e <rp'(t) 



i?i(x,t) = / <?(s)ds = e 2 <r(p(t) -p(0)) 
•/ o 

so that by (|5.23|) we get that there exists C > such that for all t G [0, T], 

sup < Cep(t). 

fe, S )G0x[O,T] 

Hence by standard parabolic estimates, there exists C > such that for all 
(x,t) G 0t, the solution W\ of Problem (A\) satisfies 

|^i(x,t)| + \VWi(x,t)\ < Cep(t). (5.25) 

Problem (A%) 

Note that by the standard parabolic estimates there exists a constant C > 
such that By definition of hi , using again (|5.23|) , we obtain that there exists 
C > such that for all (s,i) G [0,T] 2 

h 2 (y,s) <2Ce~^ do/£ - p(s)) 

<2Ce~ Xdo/2e 
AC 

-XdVe e (5-26) 
4C 

~ Xd e(K - lf Pis) 
<Ciep{s) < C'epit). 
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Thus by standard parabolic estimates, we obtain that for all (x, t) G 0? 

\w 2 (x,t)\ + \Vw 2 (x,t)\ < C'ep(t), 
which implies that there exists C > such that for all (x,t) E 0t 

\W 2 (x,t)\ + |VW 2 (x,t)| < Cep(t), (5.27) 

where we define W 2 (x, t) = Jq w 2 (x, s)ds. 
Problem (^3) 

Note that /i3(y,s) is supported in {|d(y,s)| < do}. Moreover by linearity 
we may suppose that the function /13 satisfies one of the three following 
assumptions: 

(Hi) \h3(y,s)\ < X{\d(y,s)\<ep(s)} 

(H±) \h 3 (y,s)\ < e -Md(y,s)/e± P (s)\ 

Then under respectively assumptions (Hi), (H^), we define a function h on 
R x [0, T], respectively by 

= { ^ik1r>$ft| (5-28) 

Note that \h^(y, s)\ < h(d(y, s), s), and under either of the assumptions (Hi) 
or (H2), there exists a constant C > such that for all (s,t) S [0,T] 2 



/do 
h(r,s)dr <Cep(t). (5.29) 
-dn 



edo 

-do 

Let y>(x, t) = e t w^(x, t), then in view of Problem (A3), the function cf> satisfies 



{(ft — aA(p = f in 0t 
^ = ond0x(O,T) (5 - 30) 

ov 

where f(x,t) = e t hs(x,t) and ip(x,0) = ws(x,0) = for all x G 0. We 
establish now that there exist a constant C > such that 

V(x,t)eQ T , 0<p(x,t)<Cep(t). (5.31) 

As in [2J, the solution (p(x,t) of Problem (|5.3U|) can be expressed as 

<p(x,t)= / G(x,y,t - s)f(y,s)dyds, 

J0 J\d(y,s)\<do 

with G(x, y, t) being the Green function associated to the Neumann bound- 
ary value problem in for the parabolic operator (ft — aAip. Thus for all 
(x, t) G Q,t, 

0<<p(x,t)<[ I G(x,y,t - s)e s h(d(y,s),s)dyds (5.32) 

JO J\d(y,s)\<d 
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Next we recall the following important property of G which is established 
i» 0- 

Lemma 7.6, |2j: Let V be a closed hypersurface in f2 and denote by d(x) 
the signed distance function associated with T. Then there exists constants 
C, do > such that for any function n(r) > on R, it holds that 

/ G(x, y, t)n(d(y))dy < -= / n(r)dr forO<t<T 

J\d\<d Vt J-do 

Moreover as pointed out in [2], the above inequality is uniform with respect 
to smooth variations of T and for t G [0, T]. Applying this inequality to our 
case, we deduce that there exists C > such that for all (x,y) G O 2 and for 
all < s < t < T, 

f C f d ° 

/ G(x,y,t — s)h(d(y, s), s)dy < —== \ h(r,s)dr. (5.33) 

J\d(y,s)\<d yt — SJ-do 

In view of (|5.32p and of (|5.29|) . it follows that for all x £ and for all 

te[0,T], 

< ip(x,t) < C / G(x,y,t — s)h(d(y, s), s)dyds 

J0 J\d(y,s)\<d 
"t i rd 



< C' / / h{r,s)drds 

Jo Vt - s J -d 
rt i 

< C' / ; ep(t)ds < 2C'ep{t)VT 

Jo Vt~ s 

which yields inequality (|5.3ip . 

Coming back to W3, we deduce that for all (x,t) £ Ot, 

\w 3 (x,t)\ = \e^ip(x,t)\ < Cep(t). (5.34) 

Define Ws(x, t) = f£ ws(x, s)ds, then it follows that 

\W 3 (x,t)\<Cep(t) (5.35) 

for some C > and for all (x,t) £ £It- We show now that there exist C > 
such that for all (x,t) € 0t, 

\Vd(x,t).VW 3 {x,t)\ <Cep{t). (5.36) 

Time integration of the equation in Problem (^3) on [0, t] gives 

w 3 (x,t) — u>3(x,0) — aAWs{x,t) + W 3 (x,t) = / /i 3 (x,s)ds. 
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Since 103(2;, 0) = 0, we obtain the following elliptic problem for any t 6 [0, T], 

-aAW 3 (.,t) + W 3 (.,t) = H 3 (.,t) in0 

— — L,t) = on <90 

ov 

where for all (x,t) 6 0t, 

H 3 (x,t) = / ^(x, s)ds — 103(3;, i). 



JO 

Let us define for any t £ [0, T] the functions a(.,i) as the solution of 

f -aAa(.,i) + a(.,i) = h 3 (.,t) in 

^ ^ A n 0/74 (5-38) 

and define = §ta(x,s)ds. Define similarly B(.,t) as the solution of 

f -aAB(.,t) +B(.,t) = -w 3 (.,t) in 

9B, , ^ (5.39) 

I — ,t =0 on 00 V ; 

so that by linearity 

V(x, t) £ n T , W 3 (x, t) = A(x, t) + B(x, t). 

It follows from standard elliptic estimates in view of (|5.34p that 

\B(x,t)\ + \VB(x,t)\ < Cep(t). 

Concerning a, note that the elliptic problem appearing here is the same as 
for the chemotaxis-growth system studied in [3] and in [I], with the right- 
hand-side satisfying (|5.29|) . Therefore the results stated in Lemma 4.2 in 
PQ and in Lemma 4.10 in [3] apply and prove that there exists a constant 
C > such that for all (x, t) £ T , 

\a(x,t)\ + \Vd(x,t).Va(x,t)\ < Cep(t) 

and consequently 

\A(x,t)\ + \Vd(x,t).VA(x,t)\ < Cep(t). 

This completes the proof of (OBI) . In view of (1Q5J) . (IQ71) . (jOg) and 
(|5.36p . inequality (|5.14p is now established. 

In order to prove inequality (|5. 15|) . note that using (|5.14p we obtain that for 
all (x,t) G 0t, 

\{v-v°){x,t)\ = \v (x)e- x ti mix ' s)ds - Vo {x)e- x ti m °( x ' s)ds \ (5.40) 
< C\v {x)\\ f (m-m°)(x,s)ds\<C'ep(t), (5.41) 
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where C > is a suitable constant. 

Next we have similarly that for all (x,t) G 0r, 



| Vd(x, t) ■ V(v - v°)(x, t)\ < C\e~ x m(x,s)ds _ g -A /* «."(*,.)& i 



+C\v(x,t)Vd(x,t) ■ J Vm(x,s)ds - v°(x,t)Vd(x,t) • / Vm°(>,s)ds| 

< C'ep(t) + C\v{x,t)\\Vd(x,t) ■ [ V(m-m°)(x,s)ds\ 

Jo 

+C\v(x,t) - v°(x,t)\\Vd(x,t) ■ / Vm°(x,s)ds\, 



where C,C > are suitable constants. Using ()5.4ip . (|5. 14)) and upper 
bounds on \v\ and |Vm°|, we deduce that (|5. 15 j) is satisfied. This completes 
the proof of Lemma 5.6. 

5.4 Proof of Corollary 11.31 and Theorem 11.21 

The pointwise convergence of u £ to u° in \J o<t<T ((0 \ I\) x t) when e — > 
follows from Lemma 15.41 and from (|5.7p . Next note that w e = m £ — m° is a 
solution of Problem (|5.2ip with the right-hand-side h e satisfying 

\h e (x,t)\ < hx{x,t) + h 2 (x,t) + h 3 (x,t) 

with hi, i = 1,2,3 defined as in the proof of Lemma 15.51 This shows that 
there exists C > such that 

I^ £ Hl 1 (0 t ) - Ce. 

It follows then from standard parabolic estimates and Sobolev inequalities 
that for any a G (0, 1) there exists p G (1, +oo] and C > such that 

\\w £ \\ cl+a , 1+a / 2( ) <C||n e - u°|| L p( ) 

{ ' x (5.42) 

<CW\\ LH0T) <Cev. 

Thus for any a G (0, 1), 

lim||m £ -m°|| cl+a , (1+a)/2( T) = 0. 

The expression of v £ and Vv £ in (|2.3p and (|2.7p then shows that 

lim ||v e - f O || c i +a ,(i +C( )/2(0 T ) = 
which completes the proof of Theorem 11.21 
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5.5 Proof of Theorem II. 1L Theorem 11.41 and Corollary 11.51 



The proofs are exactly the same as the proofs of Theorem 1.3, Theorem 1.5 
and Corollary 1.6 in [1] respectively, we omit the details here. 
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